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Overview 


This  final  report  summarizes  a  joint  research  project  between  the  University 
of  Delaware  Center  for  Composite  Materials  (CCM),  the  Advanced  Com¬ 
puting  and  Simulation  Division  (ACISD)  of  the  Army  Research  Laboratory 
(ARL),  and  the  Army  High  Performance  Computing  Research  Center  (AH- 
PCRC).  CCM  is  an  established  ARO  center  in  materials  science  research. 
ACISD  is  currently  supporting  an  initiative  in  virtual  manufacturing.  Un¬ 
der  the  technical  direction  of  Dr.  Kurt  Fickie,  the  research  initiative  seeks 
to  develop  a  manufacturing  design  evaluation  capability  which  will  reduce 
the  time  required  for  pr-notyping  new  weapon  systems  utilizing  fiber-resin 
composite  materials.  Such  systems  include  the  Comanche  helicopter  and 
the  Composite  Armored  Vehicle  (CAV). 

The  project  goal  was  to  complement  CCM's  leadership  role  in  experi¬ 
mental  and  analytic  composites  research  with  advanced  numerical  modeling 
methods  and  computational  resources  available  through  the  AH:'  CRC.  The 
contract  provided  graduate  support  and  trav<a  for  the  principal.  vestig<ator 
to  visit  ARL  and  CCM  over  a  two  year  period.  A  total  of  eight  (8)  trips 
were  made  by  the  PI  to  ACISD/APG  and  CCM.  Two  research  papers  were 
submitted  for  publication  in  refereed  journals,  and  a  production  software 
system  for  design  evaluation  was  jointly  developed  with  CCM  and  ACISD. 

2  Statement  of  Problem 

One  of  the  primary  design  evaluation  problems  in  manufacturing  large  com¬ 
posite  parts  is  to  determine  whether  the  part  will  be  fully  saturated  by  resin 
under  the  design  injection  conditions.  The  flow  of  a  polymeric  fluid  through 
a  porous  network  of  fibers  is  a  key  phenomenon  in  the  resin  transfer  rn.old 
filling  (RTM)  process. 

Numerical  simulation  of  the  RTM  filling  process  is  necessary  for  a  basic 
understanding  of  the  process,  and  for  addressing  production  design  needs  for 
large,  complex  parts.  Flow  modeling  is  necessary  to  predict  process  variables 
such  as  fiber  wetting  and  progression  of  the  free  surface  flow  front.  Flow 
modeling  is  useful  in  predicting  whether  a  given  model  design  is  feasible  (i.e., 
the  part  will  wet  out)  and  for  optimization  of  the  mold  design.  Some  of  the 
properties  of  mold  design  which  can  be  optimized  during  the  design  process 
include  gate  location  (a  gate  is  a  resin-injection  site),  number  of  gates,  and 
flow  rates. 


2 


Other  primary  design  evaluation  problems  include  the  prediction  of  resid¬ 
ual  stresses  and  loading  behavior  of  composite  parts.  These  problems  were 
not  addressed  by  this  effort. 

3  Summary  of  Results 

The  primary  research  need  articulated  by  ARL  and  COM  researchers  was 
the  fast  simulation  of  isothermal  (constant-temperature)  mold  filling.  The 
problem  was  that  existing  simulations  were  extremely  time-consuming  be¬ 
cause  the  time  complexity  tended  to  grow  as  O(n^),  where  n  is  the  number 
of  nodes  in  the  finite  element  mesh.  Simulation  of  a  part  with  even  a  few 
thousand  nodes  required  many  hours.  In  response  to  this  problem,  a  new 
algorithm  for  solving  the  isothermal  filhng  problem  was  developed  and  em¬ 
bodied  in  the  Fortran  code,  ISOFIL.  The  new  algorithm  involved  a  novel 
linear  algebra  approach  to  build  the  Cholesky  factorization  of  the  pressure 
stiffness  matrix  row  by  row,  rather  than  reassembling  and  refactorizing  the 
entire  stiffness  matrix  at  each  time  step.  This  algorithm  reduced  the  time 
complexity  of  the  simulation  to  in  theory,  with  an  observed  com¬ 

plexity  closer  to  O(n^).  ISOFIL  was  initially  capable  of  solving  2000-node 
problems  in  a  matter  of  a  few  minutes  on  IRIS  4D-35  workstations  com¬ 
pared  to  many  hours  for  predecessor  codes.  More  recent  improvements  to 
the  ISOFIL  algorithm  and  the  availability  of  MIPS  R-8000  64-bit  processors 
have  reduced  the  time  to  a  few  seconds.  ISOFIL  remains  a  production  code 
for  numerical  simulation  of  isothermal  mold-filling  and  is  undergoing  contin¬ 
ued  enhancement  by  ARL  and  University  of  Minnesota  researchers.  It  can 
be  used  on  an  interactive  basis  for  problems  of  moderate  size  (<  10,000) 
nodes. 

Simulation  of  non-isothermal  (temperature- varying)  flow  and  the  resin 
curing  process  requires  several  rheological  parameters  -  constants  appear¬ 
ing  in  the  rate  equations  describing  heat  generation  and  gel  conversion.  Re¬ 
searchers  have  observed  significant  differences  in  these  parameters  depending 
on  whether  they  are  obtained  in  a  batch  or  flow  environment.  In  response 
to  this  problem,  a  graduate  student,  Robert  Duh  (Mechanical  Engineering, 
University  of  Minnesota)  developed  a  parameter  estimation  methodology  for 
non-isothermal  RTM.  Duh  devised  a  Lagrangian  model  of  resin  flow  which 
captured  the  important  heat  and  mass-transfer  properties  of  RTM  in  a  one¬ 
dimensional  flow  simulation.  He  developed  an  interface  between  the  model 
and  a  nonlinear  least-squares  routine  (MINPACK)  to  estimate  rheological 
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parameters  by  fitting  the  model  to  experimental  data  on  temperature  and 
P‘  ossure  obtained  in  RTM  molds. 

Parameter  estimation  typically  requires  many  evahiations  of  the  under¬ 
lying  flow  model,  so  there  is  a  need  for  simplified  and  fast  flow  solvers.  A 
concern  is  that  simplified  models  of  flow  and  thermal  transport  may  miss  im¬ 
portant  heat  transfer  dynamics.  For  example,  a  one-dimensional  Lagrangian 
model  neglects  diffusion  in  the  axial  (flow)  direction  on  the  assumption  that 
because  the  mold  is  thin,  heat  transfer  through  the  wall  dominates  the  re¬ 
sult  (i.e.,  diffusion  through  the  wall  is  much  greater  than  in-plane  diffusion). 
A  scaling  •  nalysis  suggested  that  in-plane  diffusion  could  be  important  for 
fast  flow,'.  To  further  test  this  a^ssumption,  a  three-dimensional  model  of 
heat  trar  ler  was  (ieveloped  and  .  ^)lved  using  the  differential-algebraic  sys¬ 
tem  solvci'  DASPKF90  on  the  AilPCRC  CM-5.  A  comparison  of  three- 
dimensional  results  with  the  one-dimensional  model  indicated  that  the  as¬ 
sumption  of  negligible  axial  heat  diffusion  was  reasonable  under  standard 
manufacturing  conditions  for  RTM. 

4  Publications  and  Reports 

Two  journal  publications  were  prepared  under  this  contract,  including, 

“A  Fast  Miimerical  Method  for  Isothermal  Resin-Transfer  Mold  Filling," 
R,  Maier,  1 .  Rohaly,  S.  Advani,  K.  Fickie,  submitted  to  IntL  J.  Numerical 
Methods  in  Engineering. 

"Experimental  Estimation  of  Process  Parameters  in  Resin  Transfer  Mold¬ 
ing,  R.  Duh.  R.  Maier,  J.  Vogel,  to  appear  in  SAMPE  Quarterly. 

These  manuscripts  exist  as  separate  ARO  Interim  Technical  Reports. 


5  Personnel  and  Degrees 

Personnel  supported  by  this  contract  included  Robert  Duh.  Graduate  Stu¬ 
dent,  Mechanical  Engineering  Department,  University  of  Minnesota  (  Duh’s 
anticipated  graduation  date  is  3/96.  His  work  has  been  partially  supported 
by  this  contract),  Robert  Maier.  Staff  Scientist,  AHPCRC  (partial  support), 
and  Y.T.  Chen.  Postdoctoral  Fellow,  Chemical  Engineering  and  Materials 
Science  Department.  University  of  Minneosta  (partial  support). 
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6  Report  of  Inventions 


This  work  did  not  result  in  any  patentable  inventions.  The  ISOFIL  software 
is  currently  distributed  and  used  by  ACISD. 
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